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ABSTRACT 

The Blazhko effect is a long term, generally irregular modulation of the light curves that occurs 
in a sizeable number of RR Lyrae stars. The physical origin of the effect has been a puzzle ever 
since its discov ery over a hundred years ago. We build here upon the recent observational and theo- 
retical work of ISzabo et al. I (|2010f ) on RRab stars who found with hydrodynamical simulations that 
the fundamental pulsation mode can get destabilized by a 9 : 2 resonant interaction with the 9th 
overtone. Alternating pulsation cycles arise, although these remain periodic, i.e. not modulated as in 
the observations. 

Here we use the amplitude equation formalism to study this nonlinear, resonant interaction between 
the two modes. We show that not only does the fundamental pulsation mode break up into a period 
two cycle through the nonlinear, resonant interaction with the overtone, but that the amplitudes are 
modulated, and that in a broad range of parameters the modulations are irregular as in the observa- 
tions. This irregular behavior is in fact chaotic and arises from a strange attractor in the dynamics. 

Subject headings: stars: oscillations (including pulsations), (stars: variables:) RR Lyrae, instabilities 



1. INTRODUCTION 

A large subclass of RR Lyrae stars undergo light curve 
modulations, most of them on the time scale of some 
60 periods, although the range extends from some tens 
to some hundreds of periods. The fundamental pulsa- 
tion cyc le itself l a sts ~ 0.5 day. The effect was discov- 
ered by iBlazhkol tfl907h over a hundred years ago, and 
a number of explanations have been proposed, such as 
closely spaced pulsation modes, a modal 1 : 2 resonance, 
an oblique rotator mo del, a nonradial modal interaction, 
and c onvective cycles (|Stotherdl2010l : IMolnar fc Kollathl 
120101) . However, none of these mechanisms is without 
fault. 

The giant step toward the explanation of the Blazhko 
effect has come from the unprecedently precise and con- 
tinuous Kepler space t elesco pe observations and their 
analysis (|Szab6 et al. I I2010D . Because we are going 
to build on this work we first present a summary of 
these findings. Fourteen RR Lyr stars in their sam- 
ple undergo Blazhko modulations. Unexpectedly, three 
of these stars also display period doublings, i.e. the 
shapes of the light curves shows cycle to cycle alterna- 
tions. The depths of these alternations change during 
the Blazhko cycle. Another recent observational find- 
ing is t hat the Blazhko cycle does not repeat regularly , 
see e.g. iChadid et all ((20101 ); iKolenberg et~aTI pOlll ): 
ISodor et al.1 (|2011T l" This behaviour poses another im- 
portant constraint for the physical explanation of the 
Blazhko effect, as do the observed variations of the mean 
physical parameters o f the stars during the Blazhko cycle 
(jJurcsik et"aHl2009alF 
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On theoretical side ISzabo et al. I (|2010l ) performed a 
systematic numerical hydrodynamical modeling survey 
of RR Lyr models. They found that over a relatively 
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broad region of astrophysical model parameters the fun- 
damental pulsation is unstable and develops into a pul- 
sation with alternating cyc les. Guided by earlier work 
(|Moskalik fe Buchlerl [1990D that had shown that half- 
integer resOTiance^ancause a bifurcation to alternating 
cycles. ISzabo et al. I (|2010D did some sleuthing work. By 
computing the Floquet stability coeffic ients (|Hartmanl 
[19921: IBuchler. Moskalik fe Kovacsl H99l of the funda- 
mental pulsation and searching for resonances of the fun- 
damental mode with successive overtones showed that it 
is the 9th overtone that is in a 9 : 2 resonance and desta- 
bilizes the fundamental pulsation cycle. (It is of course a 
coincidence that it should be the 9th overtone that is in a 
9 : 2 resonance.) This 9th overtone in these RR Lyr mod- 
els turns out to be egregious in that it is a surface mode 
(dubbed 'strange mode' when it was first encountered by 
IBuchler. Yecko fe Kollathl ((19971) and IBuchler fc Kollathl 
(20011)). The 9 : 2 resonance appears in a relatively 
narrow, winding band in an Log L — Log T diagram 
((Kollath. Molnar fc Szaf^l2l)lTI ). 

The numerical hydrodyn amical simulations 
(|Kollath. Molnar fc Szabcl [20Tlh were able to pro- 
duce alternating cycles, and the ancillary analyses 
established that the 9 : 2 resonance is the cause of 
the symmetry breaking bifurcation. However, these 
hydrodynamical simulations were unable to produce 
either regular or irregular Blazhko like modulations. 

When a surface mode plays a dynamical role, the nu- 
merical hydrodynamical simulations become particularly 
sensitive to the mixing length parameters, to the zoning 
in the outer regions as well as to the surface boundary 
condition. It is therefore extremely difficult to make com- 
pletely trustworthy and robust simulations. In fact it is 
an open question whether hydrodynamics with a time 
dependent mixing length treatment of convection is able 
to produce Blazhko modulations. For these reasons it is 
of importance to resort to an alternative, complementary 
approach as we do in this paper. 
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Fig. 1. — Temporal modulation of the two amplitudes of the irregularly modulated period two pulsation. 



2. MODELING 

In the 1980's several groups developed the amplitude 
equation formalism in an effort to understand the sys- 
tematics of stellar pulsation s throughout the HR dia- 
gram (for a review, see e.g. lBuchlerlH993l ). The valid- 
ity of this powerful formalism rests on the fact that for 
mildly nonadiabatic stars, such as the RR Lyrae and the 
classical Cepheids, one can decouple the 'fast' pulsation 
and obtain equations for the slow secular behavior of the 
amplitudes and phases, akin to the Poincare - Lindstedt 
approach in celestial mechanics. 

The amplitude equation formalism is complementary 
to numerical hydrodynamical simulations. It allows one 
to understand the mathematical structure of the dynam- 
ics of the stellar models. As such it is therefore ideally 
suited for making surveys of the possible pulsational be- 
havior of a given type of star. However, when light curves 
or radial velocity curves for specific models are desired, 
one needs to resort to numerical hydrodynamics. 

The form of these amplitude equations is generic and 
depends only on the resonances that exi st between the 
excited linea r modes of pulsation (e.g. iBuchlerl 119931 ; 
lSpiegeliri995() . Amplitude equations have a solid mathe- 
matical foundation and are als o known as 'normal forms' 
in the mathematical literature (jGuckenheimer fe Holmes] 
[1981 . 

We note in passing that the first success of the am- 
plitude equation formalism in stellar pulsation theory 
was to prove the conjecture of iSimon fc Schmidt] ffl976h 
that the Hcrtzsprung progression of the Fourier phases 
of the bump Ccpheid light curves and radial veloc- 
ity curves is indeed caused by the 2 : 1 resonance 
between the fundamental pul sation mode and the lin- 
early stable second overtone (iKlapp. Goupil fc Buchlerl 
19851 iKovacs fe Buchien I1990D. We also note that 



Szabo, Kollath fc Buchlerl (|2004[) were able to map out 
efficiently the pulsational behavior of RR Lyr stars 
throughout the whole instability strip by using a judi- 
cious combination of numerical hydrodynamical simula- 
tions and amplitude equations. 

In this paper we ignore first overtone (RRc) 
and double mode pulsations (RRd) that we ad - 
dressed in e.o.lKollath. Buchler. Szabo fc Csubrvl ((2002I ) ; 
iSzabo. Kollath fc Buchlerl (j2004f ). We therefore do not 
need to include an amplitude equation corresponding to 



the first overtone in what follows. Hereafter, we label the 
fundamental mode with 'a' and the 9th overtone with 'b'. 

The amplitude equations that are appropriate for the 
modal interaction between mode a that is in a 9 : 2 reso- 
n ance with mode b (9uj n ~ 2u i b ) have already been given 
in lMoskalik fc Buchlerl (|T990h . 



da 
~dt 
db 
~dt 



(ioj a + K a - Q a \a\ 2 - T a \b\ 2 )a + c a a* 8 b 2 
(icj b + K b - Q b \b\ 2 - T b \a\ 2 )b + c b a 9 b* 



(1) 

By introducing real amplitudes, A and B, and phases 
4> a and 4> b , by defining the relative phase 

T = 9(f) b - 2<j) a 

we can cast Eqs.[T]into a set of 3 real equations 



Q a A z - T a B A )A + C a A 8 B 2 cos(r + 5 a ) 



dA , 

-§ = ( K b~ QbB 2 - T b A 2 )B + C b A 9 Bcos{T - 5 b ) 

— = 2A - 9C a A 7 B 2 sin(r + S a ) ~ 2C b A 9 sm(T - S b ) 
at 

(2) 

The k denote the growth rates of the two modes. The 
quadratic terms Q a , Q b , T a and T b appear in the am- 
plitude equations even in the absence of any resonances. 
Their imaginary parts do not have a significant effect and 
convenience we have assumed that they are all real. The 
quantities C a , C b and their phases S a , S b describe the 
resonant coupling, and 

9 

A = us b - -uj a 

is the off-resonance parameter. 

Two types of solutions exist depending on the param- 
eters of the amplitude equations: 

1. The fixed points , (dA/dt = Q 7 dB/dt = Q,dT/dt = 0). 
Eqs. ((2]) have two possible fixed points: 
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i. A = and B = 0. This fixed point corresponds 
to the trivial, nonpulsating static model. It is unstable 
when mode a is self-excited, n a > 0, i.e. inside the fun- 
damental instability strip. 

ii. A single mode fixed point with A = Aq = \Jn a /Q a 
and Bq = 0. This fixed point, when stable, corresponds 
to steady periodic pulsations in the fundamental mode 
with constant amplitude. In some regime of parameters 
this solution is unstable to an excitation of mode b, in 
which case we have 

iii. A two-mode fixed point with A = Aq and B = Bq 
(and r = To). Note that because of the 9 : 2 reso- 
nance condition, it takes two fundamental periods for 
the oscillation to repeat. This is the fixed point that 
iKollath. Molnar fc Szabol ()2011f ) found in their numeri- 
cal hydrodynamical simulations and that gives rise to the 
period alternations. 

Generally, the equations ^ are complicated enough so 
that the values of Aq and Bq must be obtained numeri- 
cally. 

In some parameter regime this fixed point can also be 
unstable, as we describe now. 

2. The oscillatory solutions for the amplitudes 
These solutions come in two types: 

i. Purely periodic oscillations in the amplitudes and 
phases, corresponding to periodic amplitude and phase 
modulations of the stellar pulsations. 

ii. Irregular oscillations (irregular amplitude and phase 
modulations) . These oscillations arise because of the ex- 
istence of a strange attractor (chaos). It is these solutions 
that correspond to irregular Blazhko modulations. 

Depending on the values of these parameters all these 
solutions show up in the numerical integration of the 
amplitude equations. 

The amplitude equations have too many parameters 
to make an exhaustive study of their solutions possi- 
ble. However we can impose some constraints. By 
placing ourselves inside the fundamental mode instabil- 
ity strip we can assume that the growth rate K a > 0. 
The quadratic terms in Eqs. @, viz. Q a , Qb, T a and Tj, 
are generally positive. We can scale the time by setting 
n a = 1 and scale the amplitudes A and B by respectively 
setting Q a = 1 and Qb = 1. 

In the following we show an example of an irregular 
amplitude modulation. The values of the remaining pa- 
rameters are chosen to be Kb = - 0.1, T a = 10.0, Tb = 
1.0, C a = 0.75, C b = 11.0, S a = 0.5, 6 b = 1.2, A = -2. 

The single mode fixed point is located at Aq = 
1 (by construction). The Floquet stability coef- 
ficients corresponding to an excitation in mode b 
(jMoskalik fc Buchlenll990h 

T h = K b - T h Al ± ^{C 2 h Af - A2) 

have the values -12.034 and 9.8343, and the single mode 
is therefore unstable to an excitation of mode b. 

The two mode fixed point is found to be located at 
A x = 0.7979, B x = 0.1918, Ti = 0.1955, and its stability 
eigenvalues are calculated to be uq = —3.3135 and a± — 
0.2057 ± i 2.7913. This fixed point is therefore unstable 
as well and an oscillation arises in the amplitudes. 

The temporal variation of the amplitudes A{t) and 




0.85 0.9 . n 0.95 



Fig. 3. — First return map on the successive amplitude maxima 
AJJj of amplitude A, namely A^^x vs. A r ^ Hax 



B{t) is displayed in Fig. Q] with a blowup in the right 
panel. One notes that, in this example, the modula- 
tions in the amplitudes are irregular, akin to the observed 
Blazhko cycles. 

In RR Lyr stars the relative growth rates kq x Pq are 
typically of the order of a percent. The amplitude mod- 
ulations occur on time scale that is set by n a that is 
therefore much longer than the fundamental pulsation 
period. In our example in which we have scaled k to be 
unity, the scaled 'Blazhko period' comes out to be ~ 2.5 
as the figure shows, corresponding thus to ~ 250 periods, 
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in the ball park of the observed Blazhko cycle periods. 
We note though that this period can easily be shortened 
or lengthened with a finctuning of the parameters in the 
amplitude equations and one can get the whole range of 
Blazhko 'periods'. 

Admittedly we have chosen the parameters in our am- 
plitude equations in a somewhat ad hoc way, but we find 
that the range of parameters over which irregular mod- 
ulations occur is quite broad. 

3. NATURE OF THE IRREGULAR MODULATIONS 

The reader may wonder how such irregular modula- 
tions arise. Because of this spiral-node nature of the 
fixed point, the trajectories in the complex amplitude 
(a, b) space spiral away from the fixed point in the plane 
of the spiral roots (<r±), get bent around in the nonlinear 
regime to be attracted back to the fixed point along the 
direction of the real, stable root (oo). Under broad pa- 
rameter conditions this results in irregular oscillations. 
This behavior is displayed in a B(t) - A(t) phase plot in 
Fig. [5] for our example. 

At first sight one might think that the dynamics of this 
strange attractor is a of a stretch-and-fold type similar 
to the well known Rossler band. As a further test, we 
have constructed a first return map with the maxima of 
amplitude A, i.e. a plot of the successive pairs of the 
A amplitude, vs. A^ ax , shown in Fig. (fig3). It 

exhibits a tent like str ucture more akin to that of the 
Lorentz attractor (e.g. iThompson fc Stewart! fl986l ). al- 
though somewhat more complex because of the split. For 
reference we have drawn a diagonal line that intersects 
the map at two very closely spaced points. At both points 
the slope of the map is greater than unity ( — slope — ll) 
which establishes the presence of a strange attractor. 

Chaotic behavior had previously been detected 
in the puls ations of the RV Tauri cla ss star 
R Scuti (IBuchler. Serre. Kollath fc MatteJ 119951: 
iBuchler. Kollath. Serre fc Matteilll996l) . It hadalso been 
predicted to be the cause of irr egular pulsations in long 
period W Vir stellar mod e ls (IBuchler fc Kovacsj IT9871: 
Kovacs G. fc Buchler ,T.Rl 119881 : IMoskalik fc Buchlerl 
19901) . It is interesting that the irregular modulations 
in the Blazhko RR Lyr stars would also be caused by 
a chaotic dynamics. There is however a big physical 
difference between the two cases. In RR Lyr the 
amplitude modulations occur on a time scale that is 
long compared to the 'period' of the pulsation, because 
this time scale is set by the pulsational growth rates 
that are small in these non adiabatic stars. In contrast, 
the W Vir and RV Tau stars have growth rates that are 
comparable to the pulsation frequencies and therefore 
allow the modulations to occur over a cycle, resulting in 
irregular pulsations (usually called semi-regular in the 
astronomical literature). 

4. OTHER RESONANCES 

In this paper we have only considered the modal cou- 
pling between the 9th overtone and the fundamental 
mode through the 9 : 2 resonance. 

An inspection of the p e riod r atio diagram (Fig. 3 of 
IKollath. Molnar fc Szabol (|2011l )) shows that the spac- 
ings between the curves are close to 0.5, and that there- 
fore there is a plenitude of additional near resonances, 



e.g. 

UJg — LUy = LUq and UJiq — LUg, = Uq. 

Furthermore there are also several integer resonances, 
e.g. 

u>j = 4luq and u>5 = 3cjo, 

although this latter resonance may not be important be- 
cause of the strong damping of mode 05. The mode 
couplings that arise from these additional resonances ob- 
viously can cause additional features and complexity in 
the light curves. For the time being we have ignored 
them. Observations of RR Lyr stars indicate that the 
Blazhko effect can be more complicated in some RR Lyr 
stars than that explored i n this paper ( e.g. m ulti-periodic 
modulation as found by iSodor et al.l (|201lD ). In these 
cases it may necessary to include these additional reso- 
nances in the description. 

5. CONCLUSIONS 

IKollath. Molnar fc Szabol (|2011[) and iSzabo et al\ 
(|2010l) discovered with the help of numerical hydrody- 
namical simulations that the Blazhko effect is most likely 
associated with the half-integer (9 : 2) resonance between 
the fundamental pulsation mode and an overtone that 
destabilizes the fundamental RR Lyr full amplitude pul- 
sation. Because of the half integer nature of the reso- 
nance (9Pg = 2Po) h takes therefore two fundamental 
periods for the pulsation to repeat, hence the occurrence 
of alternating cycles. The hydrodynamical modeling in- 
deed found this symmetry breaking to give rise to cy- 
cle to cycle alternations, i.e. short term variations. The 
Blazhko effect, however, is a long term irregular ampli- 
tude modulation of the pulsations which hydrodyamical 
simulations have not produced so far. 

In this paper, using an entirely different approach, 
namely the amplitude equation formalism, we demon- 
strate that irregular amplitude modulations can occur 
quite naturally as a result of the nonlinear, resonant 
mode coupling between the 9th overtone and the fun- 
damental mode. The phenomenon occurs over a broad 
range of physical parameters and is therefore quite ro- 
bust. Furthermore we find that the range of 'periods' of 
the Blazhko like amplitude modulations are in concor- 
dance with the observed ones. 

It is important to emphasize that the same half integer 
resonance, responsible for the destabilization of funda- 
mental RR Lyr pulsations, is also capable of producing 
period doubling or amplitude modulations, depending on 
the coefficients in the equations. 

However, observations may have other surprises in 
store and the Blazhko effect may turn out to be more 
complicated, so that this simple 2 mode coupling may 
not account f or all its complexit y. This is the case for 
CZ Lacertae (jSodor et al.1 l)2011| )). where the Blazhko- 
modulation is multi-periodic. There exist, in effect, ad- 
ditional modal resonances or near resonances that may 
need to be added in the amplitude equation description, 
and that may then lead to more complicated light curves. 

Finally, now that the fact that period alternations 
can occur in classical variable stars such as RR Lyr 
stars has been accepted by the observational commu- 
nity, time might be ripe to also target BL Her stars 
in which theory had actual ly predicted such alterna- 
tions more than 15 years ago (Buchler fc Moskalik \fl99% 
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IMoskalik fc Buchler lll993HBuchler fc B^hl^[l99l . but 
which at the time had been received very skeptically at 
best. It is true that a recent unpublished Fourier analysis 
by Buchler and Moskalik of the OGLE data of BL Her 
stars does not reveal any conclusive evidence for alterna- 
tions, but then that was certainly also the case for the 
pre- Kepler of RR Lyr. 
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